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We discuss a system of dilute Bose gas confined in a layered structure of stacked square lattices 
(slab geometry). A derived phase diagram reveals a non-monotonic dependence of the ratio of 
tunneling to on-site repulsion on the artificial magnetic field applied to the system. The effect is 
reduced when more layers are added, which mimics a two- to quasi-three-dimensional geometry 
crossover. Furthermore, we establish a correspondence between anisotropic infinite (quasi three- 
dimensional) and isotropic finite (slab geometry) systems that share exactly the same critical values, 
which can be an important clue for choosing experimental setups that are less demanding, but still 
leading to the identical results. Finally, we show that the properties of the ideal Bose gas in a 
three-dimensional optical lattice can be closely mimicked by finite (slab) systems, when the number 
of two-dimensional layers is larger than ten for isotropic interactions or even less, when the layers 
are weakly coupled. 
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I. INTRODUCTION 



Systems of dilute bosonic gases confined in optical lat- 
tices are ideal toolboxes for testing theoretical models 
and their solutions [T]. However, investigation of their 
properties under external magnetic field seemed to be 
precluded, since the particles used in the experiments are 
uncharged atoms and thus, are not directly affected by 
the magnetic field. Fortunatelly, development of various 
experimental techniques allow to investigate ultracold 
systems that are described by exactly the same Hamil- 
tonians as the ones interacting with external magnetic 
field. Furthermore, they lead to appearance of vortices in 
Bose-Einstein condensates (BEC), which is a hallmark of 
a superfluid in a magnetic field. One of those techniques 
results from equivalence between the Lorentz force and 
the Coriolis force. As a result, rotation of the condensate 
(usually rotation of the masks with set of holes located 
in the laser beams producing quasi-two-dimensional ro- 
tating optical lattice) acts as a synthetic magnetic field 
(SMF) |2H1]- However, this approach puts limit on max- 
imum rotational velocity, thus large SMF (e.g. required 
for quantum Hall physics) cannot be reached. To over- 
come those difficulties, imprinting of the quantum me- 
chanical phase is used, which is based on superimpos- 
ing of an external potential on a BEC (e.g. by applying 
rotating magnetic field) [SJ [B]. Another approaches are 
stirring a BEC with laser beam ffl, or using additional 
Raman lasers to coherently transfer atoms from one inter- 
nal state to another. This induces a non-vanishing phase 
of particles moving along a closed path, which simulates 
magnetic fiux through the lattice |s8, Additionally, 
rectification of the magnetic field in the optical lattice by 
using a superlattice allows to ensure that each plaque- 



tte acquires the same phase, thus simulating the uniform 
magnetic field for any value of phase between and tt 

A ground state of neutral atoms in an optical trap- 
ping potential can be either superfiuid (SF) or a Mott- 
insulator (MI). The zero-temperature coupling {t/U) vs. 
chemical potential (^/C/) phase diagram contains char- 
acteristic lobes marking a quantum phase transition be- 
tween SF and MI states; t is the matrix element for tun- 
neling between adjacent lattice sites and U is the on-site 
energy cost for multiple occupancy. For average number 
of particles per site equal to one, the transition occurs 
at the ratio (i/t^)crit' which is strongly dependent on 
the geometry of the system. Recently, many very pre- 
cise calculations and measurements have been performed 
to obtain correct value of the {t/U)^^-^^. For the three- 
dimensional {iD) optical lattices, most methods converge 
to {t/U)J^^^ = 0.03 |ITJ[I2], but lower-dimensional cases 
are more demanding because of growing infiuence of the 
quantum fiuctuations. In order to precisely describe ex- 
perimental results, it is necessary to develop a theory 
that contains fully tunable lattice degrees of freedom, in- 
cludes effects of the synthetic magnetic field and is non- 
perturbative, which allow to capture essential physics of 
strongly correlated system in U/t ^ 1 regime. To this 
end, we constructed a theoretical field approach, in which 
the dimensionality along c axis can be tuned by adding 
an arbitrary number of single layers along with vari- 
able tunneling between them (thus including anisotropy 
of tunnelling between and within planes). In principle, 
we can include sixty layers as in the experiments of the 
Spielman's group [13], or just a few as in the Kriiger 
setup ^14j . where the magnetic confinement along the z- 
direction was used to localize atoms in the effectively 
TV = AkBT/mujlP ^ 2 ^ A central lattice planes (T is 
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the temperature, ojz is the frequency at the bottom of 
the lattice wells and I is the lattice period that can be 
adjusted to any value higher than half of the atomic res- 
onance wavelength). The Burger group [15; developed 
method in which by increasing the nodal planes of the 
optical lattice superposed to a 2>D potential it is possible 
to follow the transition from three-dimensional system 
to an two-dimensional array. Therefore the system be- 
comes quasi-3Z) rather than purely planar. That can be 
used in order to look in the interference between planes 
and in consequence to access to spatial coherence of such 
structure. 

The outline of the paper is as follows: in Sec. II we 
introduce the model Hamiltonian and the effects of arti- 
ficial magnetic field and briefly describe our method. In 
Sec. IV, we present our results starting with the zero- 
temperature phase diagrams and its dependence on the 
system geometry. Furthermore, we show the influence 
of the synthetic magnetic field on properties of the Bose- 
Hubbard system. Also, we find a correspondence between 
systems sharing the same critical properties but differing 
with the geometry of the opticall lattice and interactions, 
which can be a helpful clue on choosing between various 
experimental setups. Finally, we conclude in Sec. V. 



II. MODEL 

In optical lattices, two main energy scales are set by 
the hopping amplitude t (the kinetic energy of bosons 
tunneling between the lattice sites), and the on-site re- 
pulsive interaction U (resulting from repulsion of multiple 
boson occupying the same lattice site). For t U, the 
superfiuid order is well established in zero-temperature 
limit. However, for sufficiently large repulsive energy U, 
the quantum phase fiuctuations lead to suppression of the 
long-range phase coherence resulting in SF to MI transi- 
tion. The critical ratio of [t/U)^ for which this transition 
occurs depends strongly on the number of bosons intro- 
duced to the optical lattice (which in theoretical models 
is often controled by a chemical potential /i). The syn- 
thetic magnetic field B (resulting either from rotation of 
the system, phase imprinting, or external electric field) 

— f' A. dl 

introduces the Peierls phase factor e j , where 
B = V X A(r), and $o = hc/e is the flux quantum, 
with A (r) being the vector potential (which can be re- 
alized experimentally, see Ref. ^), and h, c and e - 
Planck constant, speed of light and charge of electron, 
respectively. Thus, the system can be described by the 
following quantum Bose-Hubbard Hamiltonian [TBI [T7] 

n = ^ y^n^ {rir - 1) 

r 

— ^ t^^'e ■'^^J '''aj.a^' — ^ n^, (1) 

(r,r') 



where aj. and a^i are for the bosonic creation and anni- 
hilation operators that obey canonical commutation re- 
lations [a^,a^,] — Srr'; n-r ~ al.ar is the boson num- 
ber operator on the site r. Here, (r, r') denotes sum- 
mation over the nearest-neighbor sites. Furthermore, 
tj-r' is the hopping matrix element with the dispersion 
tk = 2t|| ^cos fcj; + cos /cj, + cos fc^^ where t± is the 

hopping between layers and t|| within the planes. Since, 
we are interested in investigating the influence of the 
lattice geometry on the system properties, we consider 
a stack of an arbitrary number (L) of two-dimensional 
planes coupled with t±. As a result, the values of and 
ky are continuous {k^^y = — tt, . . . , tt), while k^ is descrete 
(A:^ = where I = 0, . . . , L — 1). Also, we allow for c- 
axis anisotropy, which is a ratio of inter-plane to in-plane 
hopping ri = t±/tii. 

A boson hopping around a lattice cell of the area of 
A will gain an additional phase 27r/ resulting from the 
synthetic magnetic fleld, where / = ABe/2'iTh. As a 
result, the properties of the system will be periodic with 
a period corresponding to 1//. As a result, the periodic 
potential leads to splitting of Landau levels into integer 
number q of sub-bands. Of special interest are the values 
of the SMF which correspond to rational numbers of / = 
p/q = 1/2,1/3,1/4,... {p is an integer), since for those 
values the energy spectra and the density of states can be 
obtained exactly, although it is analytically feasible only 
for small values of g. In the present paper, we present 
new results for / = 1/8 and / = 3/8, which up to now 
have been analitically unaccessible (see, Ref. fT5]). Since, 
all properties of the Hamiltonian Eq. ([T|) are invariant 
under f ^ —f and also under / ^ / + 1, it is sufficient 
to consider / in the range < / < 1/2. In solid state 
physics obtaining / = 1/2 in experimental setup would 
require magnetic field of the order of lO^T. In optical 
lattices however, due to larger lattice spacings (and larger 
A consequently) , much smaller values of B of the order 
of lO^^T are required. This makes investigation of the 
above-mentioned rational values of / reasonable. 

To proceed, we rely on the quantum rotors approach. 
The method is extensively described in Refs. f l21 ITfl \W\ 
|20j, so here we only summarize its main points. We use 
the functional integral representation of the model with 
bosonic operators becoming complex fields (r) (where 
r is imaginary Matsubara's time). The most important 
element of our method is a local gauge transformation to 
the new bosonic variables: 

ai (t) = 6j (t) exp (r)] . (2) 

This allows to cast the strongly correlated bosonic prob- 
lem into a system of weakly interacting bosons, sub- 
merged into the bath of strongly fluctuating gauge po- 
tentials on the high energy scale set by U. It also allows 
to formulate the problem in the phase representation, 
which is the best suited to describe MI to SF transition, 
since it is governed by phase fluctuations. As a result, 
the superfluid order parameter can be written as: 
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(3) 



where non-zero value of '0B=(exp[i0i (''")]) results from 
phase ordering and bo is the amplitude of the bosonic 
field: 



2 1- 



The coeflicient 4 + 2(1 
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) ^ is an effective number of 
nearest neighbors averaged over all lattice sites. In the 
zero-temperature limit we arrive at the equation for the 
phase order parameter: 
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with 



(5) 
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where frac (x) = x — [x] is the fractional part of the num- 
ber and [x] is the floor function which gives the greatest 
integer less then or equal to x and the phase stiffness 

Since the number of layers L in the system is finite, 
the summation in Eq. [5] runs over discrete values of 
and continuous values of k,j; and ky. However, because 
density of states of a single layer under SMF {pf ) is 
known, we explicitly derive the density of stated of the 
whole stack of L coupled planes (for calculation details, 
see Ref. [2T]): 



r/ cos kz 



(7) 



As a result, the critical line equation (-05 = 0) including 
the effects of SMF and c-axis anisotropy reads: 



1 = 



1 



(8) 



with ^0 being the half-width of the band dispersion for 
selected value oi f = p/q. 



III. RESULTS 

A. Phase diagram 

The Eq. ([s]) allows us to calculate the zero- 
temperature phase diagram of the investigated Bose- 
Hubbard model from Eq. ([T]) as a dependence of critical 
interaction on the chemical potential, SMF, number of 
layers and c-axis anisotropy: 



(9) 




Figure 1: (Color online) The zero-temperature phase diagram 
of a stack of square lattice planes (number of particles per 
lattice site is ns = 1 inside the first and ns = 2 and 3 
in second and third lobe, respectively) with magnetic field 
/ = 1/2 for various number L of layers and 77 = 1. Within 
the MI phase the phase order parameter ~ 0- 



The diagram is plotted in Fig. [T]for different number of 
layers L, synthetic magnetic field / = 1/2, in isotropic 
case (?7 = 1). In the weak coupling limit (ty ^ U), the 
kinetic energy dominates and the ground state is a de- 
localized superfluid, described by nonzero value of the 
superfluid order parameter '^b ^ 0- On the other hand, 
in the strong coupling regime (iy ^ U) the phase fluc- 
tuation become signiflcant and the long-range order is 
destroyed leading to a series of MI lobes with flxed in- 
teger filling ub = 1, 2, ...jT21 IS]. A single-layer system 
(L = 1) has a simple square (two-dimensional) geom- 
metry, which results in the phase diagram with charac- 
teristic narrow-edged lobes. As the number of layers is 
being increased, the tops of the lobes become smooth 
and their maxima deviate towards lower values of the 
chemical potential /i (this effect is also clearly presented 
in Ref. ^T). As a result, the phase diagram becomes 
similar to the one of a cubic (three-dimensional) system. 
It is important to notice, that also in the presence of 
the synthetic magnetic field, the phase diagrams of the 
finite L system becomes indistinguishable from the in- 
finite (cubic) one for L as small as 10. For example, 
for ten layers, the value of x\J^ is 102.2% of the infinite 

system z^^g and for sixty - x^^ is close to 100.1% of 
x^jyg. The phase diagram obtained from the quantum ro- 
tor approach (QRA) was verified several times [17J [T5] 
by comparison with the very precise numerical methods: 
quantum Monte-Carlo |11| and diagrammatic perturba- 
tion theory [22_. The QRA phase diagrams were also 
compared to an analytical works: mean-field theory [2S\ 
and Fade analysis [531 . From the above, it is expected 
that QRA accurately captures the low energy physics of 
the Bose-Hubbard model in the more demanding case 
where the system is placed under an artificial magnetic 
field. 
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Figure 2: (Color online) Ratio of critical coupling x^/x^ of 
the tip of the first lobe [see, Eq. of the system without 
(/ = 0) and with synthetic magnetic field / as a function 
of number layers L for isotropic (top) — 1 and anisotropic 
(bottom) 77 = 0.1 case, respectively. 



B. Effect of the synthetic magnetic field and the 
system geometry 

The transition, experimentally seen in the time-of- 
flight images (the presence of sharp peaks has been con- 
sidered as an unequivocal signature of superfluidity in the 
Bose system), occurs rather rapidly with increasing lat- 
tice depth. Because the experimental parameter Vq/Eu 
(Vq is the maximum value of the lattice depth), depends 
logarithmically on U/t, the small changes of the dimen- 
sionless depth of the optical lattice can cover a wide range 
of the phase diagram. The phase coherent Bose gas can 
be also driven into the Mott insulating phase by apply- 
ing the synthetic magnetic field. The effect of the SMF is 
presented in Fig. [2] The long range order is suppressed 
by the phase changes imposed on the bosonic wave func- 
tion and this suppression has a non-monotonic character 
strongly depending on the topology of the system. The 
Mott insulating phase becomes more stable, which is as 
expected since the magnetic field should localize parti- 
cles. In the single-layer system, the effect of the SMF 
is the most pronunced and this decreases with growing 
number of layers L. By adding more layers the global 
coherence of the system is restored, because growing di- 
mensionality entails the suppression of quantum fiuctua- 
tions effects. Here, the convergence of properties of the 
finite system to those of the infinite (cubic) one is much 
slower, although also non-trivially dependent on / (for 
some values like / = 1/8, 1/6 and 3/8 seems to be much 
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Figure 3: The evolution of the Hostadter butterfly with in- 
creasing number of layers L in the isotropic case of tunneling 
ratio t± — Adding more layers simply results in additional 
fractal patterns and very complex structures of L butterflies 
emerging. 



more pronunced than for / = 1/4,1/3 and 1/2). When 
the system is anisotropic (see, the bottom plot in Fig. 
U, the convergence is much faster, but still depending 
on the specific value of /. 

In the non-interacting ([/ = 0) Bose-Hubbard model 
under the synthetic magnetic field, the energy spectrum 
is known as the Hofstadter butterfly |25] (see, the top 
plot for i = 1 in Fig. [3]). The band width is strongly 
dependent on / and exhibits a self-similar gap structure 
for rational values of /. It has been observed, that in 
the mean-fleld approach, the critical hopping xj- roughly 
follows the bandwidth of the Hofstadter's butterfly [23 . 
However, in our results the non-monotonicity of Xf is not 
exactly following the Hofstadter butterfly bandwidth, es- 
pecially for lower values of /. Also, we calculate the band 
for multiple-layer system and present it in Fig. [3] It is 
clear, that increasing the number of layers (for L > 1) 
does not influence the width of the band, while the Xq /x^ 
in Fig. [2]evidently changes while converging to the three- 
dimensional (cubic) case. We would like to note that, 
both densities of states for square system in the mag- 
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Figure 4: (Color online) The evolution of normalized critical 
ratio Xq jx^ with number of layers for magnetic field f = 1/2 
and / = 1/4 with different values of the tunneling ratio r] = 
t±/t\\. Dashed lines show correspondence between systems 
with finite and infinite L. 



netic field pf and DOS of L-layer stack {ViO £^re ob- 
tained exactly in analytical form with any approxima- 
tions. The Hofstadter's band structure has also been in- 
vestigated in context of quantum graphs and the integer 
quantum Hall effect [26J, where similar energy spectra for 
two-dimensional to three-dimensional system crossover 
have been observed (non-symmetricity of the spectra re- 
sults from their representation as functions of wave vec- 
tors instead of energy) . 



C. Correspondence between finite L-layer and 
infinite systems 

Although, conventional optical lattices typically pos- 
sess a uniform tunneling matrix elements, by changing 
lasers intensity along z axis one can control the tunnel- 
ing between adjacent layers. Sometimes the tunneling 
rate between adjacent sites is just negligible on the time 
scale of experiments [14]. On the other hand, using a 
mask in the novel holographic methods [27, one can also 
reproduce almost arbitrary potential. Such fully control- 
lable optical environment allows us to restrict the spatial 
(along z axis) degree of freedom of the particles and in- 
troduce c-axis anisotropy to the system. In consequence, 
the tunneling between planes can be different from that 
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Figure 5: (Color online) The evolution of normalized critical 
parameter with number of layers without magnetic field 
and for / = 1/4 with different values of the tunneling ra- 
tio rj — t^/t\\. Dashed lines show correspondence between 
systems with finite and infinite L. 



in the single layer tj^ 7^ i|| or even completely suppressed. 
However, since some experimental setups can be more 
convenient to use than others, once can try to establish a 
correspondence between number of layers and anisotropy, 
i.e. find values of L and rj that result in the systems that 
share the same properties in time-of-fiight experiments, 
thus are interchangable. The results are presented in 
Figs. |4] and [5] for two setups: when the effect of the field 
is compared to the system in zero field {xq/xj, see Fig. 
and when the bare critical value of Xj is monitored 

L 



14 

(see 



Fig. [5|. In both cases, we calculate Xq/x^ or 



as a function of finite number of layers L for isotropic 
(rj = 1) system. Then, we determine the anisotropy ratio 
r] for the infinite cubic system (L — >■ 00), which results in 
the same x^/x^ or x^ values. For example, from Fig. |4] 
it is clear that 4-layer (quasi 2D) isotropic system corre- 
sponds to 3D cubic system with anisotropy ratio equal to 
■q = 0.755 (/ = 1/2, top plot). If the field value / = 1/4, 
a 3D system of similar anisotropy (77 — 0.708) shares the 
properties with 10-layer isotropic one. Fig. [sjshows, that 
the c-axis anisotropy is less important when the magnetic 
field is turned off (values of rj are closer to isotropic rj = 1 
case) . 
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IV. CONCLUSIONS 

We have calculated critical properties of the layered 

structure of strongly interacting bosons under a synthetic 
magnetic field. Derived phase diagram reveals a non- 
monotonic dependence of the ratio of tunneling to on-site 
repulsive interaction on artificial magnetic field applied 
to the system. By calculating analytically the density 
of states (cases with / = 1/8 and / = 3/8 have not 
been presented in the literature) for several values of the 
magnetic field we were able to accurately predict the evo- 
lution of the system towards the Mott phase. The effect 
is reduced when more layers are being added, i.e. during 
the two- to quasi three-dimensional geometry crossover. 



Furthermore, we have established a correspondence be- 
tween anisotropic infinite (quasi three-dimensional) and 
isotropic finite (slab geometry) systems that share ex- 
actly the same critical values, which can be an impor- 
tant clue for choosing experimental setups that are less 
demanding, but still leading to the identical results. Fi- 
nally, we have shown that the properties of the ideal Bose 
gas in three-dimensional optical lattice can be closely 
mimicked by finite (slab) systems, when the number of 
two-dimensional layers is larger than ten or even less, 
when the layers are weakly coupled. 
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